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ABSTRACT 

We explore the nature of ‘hills’ observed on the solar surface which had previ- 
ously been attributed to Rossby waves. We investigate the solar hills phenomenon 
by analyzing the output from a synthetic model based solely on the observed solar 
photospheric convection spectrum. We show that the characteristics of these hills 
can be explained by the corrugation of the surface produced by the radial flows 
of the convection. The hills in our simulations are dominated by supergranules, a 
well-known component of solar convection. Rossby waves have been predicted to 
exist within the Sun and may play an important role in the dynamics of the so- 
lar interior, including the Sun’s differential rotation and magnetic dynamo. Our 
study suggests, however, that the hills observed at the solar limb do not confirm 
the existence of solar Rossby waves. 

Subject headings: convection — Sun: helioseismology — Sun: photosphere — 
Sun: rotation — turbulence — waves 


1. Introduction 

It is well known that the shape of the Sun is slightly non-spherical, its oblate shape 
exhibiting bulging at the equator and flattening at the poles due to rotation. Observations 
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using the MDI and RHESSI instruments both provide up to date evidence of this oblateness. 
More detailed studies of the collected data by Kuhn et al. (1998, 2000) reveal further varia- 
tions in the shape of the solar limb, with evidence for small scale ‘hills’ on the surface. These 
hills have been characterised as being 100 m high and uniformly spaced by nearly 90,000 km 
over the solar surface and were thought to be a manifestation of solar Rossby waves (Kuhn 
ct al. 2000). 

Rossby waves (inertial waves whose restoring force is the latitude-dependent strength 
of the Coriolis force) were first proposed by Gilman (1969a, b) to operate in the interior 
of the Sun. He suggested that these waves contribute to a magnetohydrodynamic dynamo 
that produces solar magnetic fields, including their periodic polarity reversals. Rossby waves 
could also be important in maintaining the structure of the Sun’s differential rotation (see 
Kato k Nakagawa 1969; Tikhomolov 1998; Wolff 1998). These theoretical studies suggest 
that Rossby waves may be connected to a large array of phenomena such as the solar surface 
structure, the Sun’s magnetic cycle, the initiation and propagation of magnetic waves, and 
the dynamics of the solar atmospheric environment responsible for “space weather” . It has 
been argued that similar effects also exist in other stars (Papaloizou k Pringle 1978) and 
gaseous giant planets (Schaeffer k Cardin 2005). Rossby waves were previously observed in 
the Earth’s oceans (Chelton k Schlax 1996; Cipollini ct al. 2001) and atmosphere (Kursinski 
et al. 1996) as low-amplitude long wavelength features. 

In their discovery of the Rossby wave ‘hills,’ Kuhn et al. (2000) argued against the possi- 
bility that these hills were produced by supergranulation. Supergranulation is a photospheric 
convective phenomenon, consisting of cells typically 30,000 km in diameter that last for about 
a day (Leighton et al. 1962). Supergranulation cells are distributed fairly uniformly over the 
solar disk. Their boundaries are where much of the magnetic flux through the photosphere is 
concentrated — forming the magnetic foundation of the chromospheric network. On global 
scales, the transport of magnetic flux across the solar surface by supergranular flows plays 
a critical role in the evolution of the large-scale surface magnetic field. Such surface field 
patterns ultimately determine the photospheric footpoints of magnetic field lines as well as 
the chromospheric and coronal field topology, which is key for their lion-radiative heating 
and dynamics. It has been found that the radial flow in typical supergranules is about 30 
m s -1 , corresponding to only 10% of their associated horizontal flows (Hathaway et al. 2000, 
2002). These radial flows should produce a series of bumps on the surface of the Sun but 
Kuhn et al. (2000) argue that the spacing would b’e too close and that the random forcing 
of these flows would not produce the features they see. 

In this Letter, we will examine the nature of hills on the solar surface previously at- 
tributed to Rossby waves (Kuhn et al. 2000), which will lead us to the conclusion that their 
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existence can indeed be explained by the presence of supergranules. In §2, we describe the 
employed methods and results. Our conclusions are given in §3. 


2. Methods and Results 

We take the following approach. Photospheric velocity data are produced using the 
method described by Hathaway (1988) that is based on the construction of synthetic data 
which realistically match observed photospheric velocity fields (Fig. 1). These simulations 
have successfully been used in numerous previous studies (Hathaway et al. 2000, 2002, 2006). 
The radial component of the velocity field produces variations in photospheric heights which 
are then analyzed in a manner equivalent to that employed by Kuhn et al. (2000) in their 
analysis of the solar hills. 

The limb signal data are produced from simulated patterns of the radial velocities on 
the solar disk following the method of Hathaway (1988). The radial velocities are generated 
by spectral coefficients R( m and the spherical harmonics W m ($, 4>). The magnitudes of the 
spectral coefficients are obtained by matching the Doppler velocity patterns observed from 
the MDI instrument on the ESA/NASA SOHO Mission with the Doppler velocity patterns 
produced in the simulations (Hathaway et al. 2000, 2002). A simplified photospheric con- 
vection spectrum can match these data over a wavenumber range of 1^1 to 1000 using two 
Lorentzian profiles, one for supergranules with a peak at £ ~ 110 and another for granules 
with a peak at £ ~ 4000. The parameters used to construct these component spectral pro- 
files are adjusted to attain a close match between the simulated and observed photospheric 
convection spectra (Fig. 2). 

The spectral coefficients for the radial velocity are converted to spectral coefficients 
for height by dividing the velocity amplitude by the square root of the wavenumber. This 
relation is adapted from the theory for up-flowing gas forming a ‘forced plume’ subject to 
buoyancy with mass flux and momentum conservation (Turner 1973). Realistic numerical 
models for solar granulation (Stein & Nordlund 1998) indicate that these smaller convective 
features have diameters of about 1000 km, radial velocities of about 3000 m s -1 , and produce 
surface height variations of about 30 km. Using the characteristic size (30,000 km) and radial 
velocities (30 in s _1 ) for supergranules indicates that supergranules should produce ‘hills’ 
with heights of about 2 km - more than an order of magnitude larger than the hills reported 
by Kuhn et al. (2000)! 

We produce height maps of 4096 by 1024 elements (longitude by latitude), representing 
the full solar surface in a similar manner to how the Doppler velocity maps had been produced 
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in previous studies (Hathaway 1988; Hathaway et al. 2000, 2002. 2006). Two latitudinal pixel 
columns separated by 2048 longitude pixels (i.e. meridional columns separated by 180° in 
longitude) are extracted from the height map array. These columns are reallocated into 
a single one-dimensional array representing the limb height all the way around the Sun 
starting from the west (as seen in the sky) and proceeding counter-clockwise through north. 
The center of the right-hand column represents the western limb, and the data north of that 
represent the northwestern limb (ending at the north pole). The data from the left-hand 
column are joined on (north to south) representing the eastern limb then the remainder of 
the first column is joined onto that (south pole to western limb). This is repeated for each 
time step (a solar rotation/4096) by stepping through the entire height map one longitude 
step at a time. (Note that halfway through the extraction, the western limb data column 
wraps around back to the far left of the map.) These limb height strips are combined into a 
two-dimensional array, so they represent a series of limb height data around the solar limb for 
each 9 minute time interval, or conversely, a series of limb height data at particular positions 
around the limb over one whole 27 day rotation sampled at a cadence of 9 minutes. A set 
of thirty such 27 day time-series are produced, analyzed, and averaged for the purposes of 
noise reduction. 

Additional data sets are produced representing limb heights over 34 days for evolving 
supergranule patterns. The patterns are evolved by adding random changes in phase to the 
spectral coefficients to give a supergranule lifetime of around 24 hours, thereby removing 
any long-range coherence in the time-series. These data are then reformatted into 34 day 
time-series with a 12 minute cadence (matching the cadence used by Kuhn et al. 2000). 

Each limb scan strip is convolved with a 4 pixel wide filter and the data sampled every 
fourth pixel, to create 512 limb position bins similar to those produced by Kuhn et al. (2000). 
Taking the FFT of each of these 512 pixel spatial strips produces a spatial spectrum for each 
time sample. Subsequent sampling along the time strips and taking the FFT of each strip 
similarly produces a temporal spectrum for each spatial frequency. This results in two- 
dimensional arrays which are the 2-D FFT of the original spatio-temporal data. Multiplying 
these 2-D FFT arrays by their complex conjugate gives power spectra for the limb height over 
a range of spatial and temporal frequencies. The power spectra are summed over all spatial 
frequencies (£ — 1 to 255) and then multiplied by the temporal frequency to produce the 
spectra shown in Fig. 3. It is clearly seen that an excess of power is present at a frequency 
of around 40 pHz (slightly higher than the frequency of the peak seen by Kuhn et al. 2000). 
This feature is produced by the supergranules in our simulation. The smaller granules are 
not well resolved in SOHO/MDI data as the data are further smoothed and re-sampled to 
get data at 512 limb positions. 
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The simulated data can also be used to produce temporal spectra of the limb signal as 
a function of limb position angle using a process similar to that of Kuhn et al. (2000). The 
data are smoothed in the spatial direction with a 256 bin running mean and in the temporal 
direction with an 8.3 day running mean, and the residual between the smoothed and original 
data calculated. The temporal FFT of this residual power spectrum is taken at each limb 
position and smoothed with a 4 pixel bin running mean for the 27 (or 34) day data stream 
(corresponding to the 88 pixel bin running mean over the 800 day sequence used by Kuhn 
et al. 2000). Multiplying this result with its complex conjugate produces temporal power 
spectra which can be plotted against position angle, as shown in Fig. 4. A ‘sawtooth’ pattern 
is seen, similar to that produced by Kuhn et al. (2000), with the maximum and minimum 
frequencies occurring at the equator and the poles, respectively. 

The temporal power spectra at each limb position are stretched in proportion to the lin- 
ear speed of rotation at the corresponding latitudes so that the maximum for each limb angle 
now lies along the frequency of maximum power at the solar equator instead of distributed 
along the sawtooth pattern. The power spectra are then averaged over all limb angles and 
plotted against the temporal frequency (Fig. 5). The non-evolving supergranule patterns 
give a peak near 30 //Hz while the evolving supergranule patterns have reduced power with 
a peak near 15 ^Hz - almost identical to the frequency obtained by Kuhn et al. (2000). 


3. Conclusions 

In this Letter, we explored the nature of ‘hills’ observed on the solar surface which had 
previously been attributed to Rossby waves. The solar hills phenomenon is investigated by 
analyzing the output from a synthetic model based solely on the observed solar photospheric 
convection spectrum (Hathaway et al. 2000, 2002). We show that the characteristics of these 
hills can be explained by the corrugation of the surface produced by the radial flows of 
the convection. The hills in our simulations are dominated by supergranules, a well-known 
component of solar convection. This result is in obvious disagreement with the conclusion 
by Kuhn et al. (2000) that the analysis of the solar limb observations described in Kuhn 
et al. (1998) constitute the detection of 100 m high photospheric ‘hills’ with 90,000 km 
separation due to the surface manifestation of Rossby waves. However, we show how the 
signatures that led to these conclusions can be created by a solar convection model based 
on a convection spectrum containing only supergranule and granule components. No wave 
phenomena whatsoever are included. Moreover, it is expected that the observational data 
should be dominated by this (larger) convection signal and we conclude that this indeed is 
the source of the observed signals. 
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We argue that what have been detected are ‘hills’ due to convection. Our model was 
also run for three scenarios, both convection components, as described above, granulation 
only and supergranulation only. On removing the supergranulation, it is found that despite 
the spectral overlap of the supergranule and granule components, the excess power features 
as seen in Figs. 3 and 5 disappear, thus granulation cannot be responsible for the observed 
signal as expected due to the small size of these features. However, the feature reappears 
for the spectrum containing only the supergranulation component. Thus we conclude that 
the signal previously attributed to Rossbv waves is most likely due to supergranulation. On 
comparing our results to the observations of Kuhn et al. (2000), however, we find that the 
latter have provided an independent technique for identifying supergranules, in addition to 
Doppler velocity maps (Hathaway et al. 2000, 2002, 2006), which may harbor future research 
potential. 

In a recent paper, Kaladze & Wu (2006) used the results of Kuhn et al. (2000) to 
explain the hill phenomena in the context of Rossby waves. As these observations can now 
be explained by supergranulation, it can be concluded that Kaladze &; Wu attempted to 
base their model on what turns out to be supergranule features. It is conceivable that what 
their model describes are indeed underlying features not yet detected on the Sun, but at 
most this work is deemed speculative and not backed up by unequivocal observational data. 

Direct observational evidence of solar Rossby waves are still yet to be confirmed. While 
our work presented in this Letter concludes that the signals reported by Kuhn et al. (2000) 
are explained by supergranulation, there may yet be Rossby wave signatures still hidden 
within the observational data reported by Kuhn et al. (1998). Data collected by the RHESSI 
instrument (Lin et al. 2002) as reported by Fivian et al. (2005) may also contain Rossby wave 
signatures. Both sets of data, however, would require the necessary filtering, including the 
removal of the supergranulation signal, to extract any signals pertaining to Rossby waves. 


This work was supported by NASA’s Office of Space Science through a grant from its 
Solar and Heliospheric Supporting Research and Technology Program. 
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Fig. 2. — Comparison between the observed photospheric convection spectrum (solid line) 
and that for the simulation (dotted line) created from granule and supergranule components. 
The spectra match exceedingly well for wavenumbers £ from 10 to 1000. The differences at 
the lowest wavenumbers are attributed to uncertainties in the instrumental artefacts found 
in the SOHO/MDI data. The peak at wavenumbers near 110 is due to supergranulation. 
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Fig. 3. — The power spectrum created by taking the spatial and temporal FFTs of the 
limb height time-series summed over wavenumber t and rescaled to similarly represent the 
plot constructed by Kuhn et al. (2000). The solid line is the average from three evolving 
supergranule simulations. The dotted line is the average from 30 non-evolving supergranule 
simulations. The excess of signal seen as a bump at a frequency of around 40 //Hz is in good 
agreement with, but somewhat higher than, that seen by Kuhn et al. (2000). The smooth 
line represents the average spectrum from the analysis of the 30 individual height maps. 
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Fig. 4. — The average power spectrum from the three evolving flow pattern data for a 
range of temporal frequencies at each angular position around the limb. The signal strength 
increases logarithmically as color changes from red to green to blue. The result is a near 
sawtooth distribution of excess power, similar to that found by Kuhn et al. (2000). This 
illustrates that by distributing similarly sized features over a rotating sphere, features at the 
equator are more numerous and pass over the limb more frequently than at the poles. 
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Fig. 5. — The spectral data from which Fig. 4 is produced are stretched for each limb 
position in proportion to its linear rotation velocity relative to that of the equator. Those 
data are then averaged over each temporal frequency bin and smoothed. The results for the 
30 non-evolving supergranule simulations display a peak near 30 juHz as shown by the dotted 
line. The results for the three evolving supergranule simulations display a peak near 15 /uHz 
as shown by the solid line. This spectral peak is similar to that found by Kuhn et al. (2000) 
from SOHO/MDI data. 




